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Abstract 

We present a detailed analysis of the results of our numerical 
study of the crystal and electronic structure of the room temperature 
organometallic ferrimagnet of general composition V(TCNE)a; with 
a: ~ 2. The results of the LSDA+C/ study show that the experimen- 
tally determined structure complies with the magnetic measurements 
and thus can serve as a prototype structure for the entire family of 
the M(TCNE)2 organometallic magnets. The results of the numerical 
study and of the magnetic experiments are interpreted using model 
Hamiltonians proposed here. This allowed us to obtain estimates of 
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the critical temperature in three- and two-dimensional regimes and to 
give an explanation of the differences in behavior of probably isostruc- 
tural V(TCNE)2 and Fe(TCNE)2 species. 

1 Introduction 



The most spectacular room temperature organometallic magnet of composi- 
tion V(TCNE)3; • ysolvent (where TCNE - 1 - stands for tctracyanoethylene 
- a well known organic electron acceptor; x ^ 2 and y depends on the type of 
the solvent) attracts a lot of attention since the time it had been synthesized 
yet in the beginning of the 1990's [1]. It is an amorphous moisture sensitive 
precipitate with the outstanding critical temperature of the transition to the 
magnetically ordered state estimated to be of ca. 400 K.^ It is higher than 
the decomposition temperature (ca. 350 K ) which singles it out among 
its numerous analogs with a variety of involved organic acceptors (Ref. [2] - 
tetracyanopyrazine - 2; Ref. [3] - 7,7,8,8-tetracyano-p-quinodimethane - 3; 
Ref. [4] - tetracyanobenzene - 4) and of the metals (Ref. [5] - iron), syn- 
thesized in the following years, since none of them manifested as fascinating 
magnetic properties as the very first V(TCNE)2 compound (see Table 1 [6]). 

Generally one has to say that not only the critical temperature, but also 
other properties of the compounds of the considered class are sensitive to 



^The magnetic momenta are spontanelusly predominantly aligned in one direction be- 
low the critical temperature. 
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the details of the preparation procedure and/or the solvent employed. For 
example, the V-TCNE compound is known in two forms. The original of 
Ref. [1] coming from the reaction of V(C6H6)2 with TCNE in the CH2CI2 
solution or obtained from V(C0)6 with use of the CVD technique exhibits 
the saturation magnetization at the zero temperature which corresponds to 
approximately one unpaired electron per formula unit and in fact is somewhat 
lower than this value. However in Ref. [7] (see also a review [6]) a V(TCNE)2 
compound was reported containing tetrahydrofurane as a solvent with the 
saturation magnetization almost twice as strong compared to that of the 
original compound. For that reason hereinafter we shall refer to these two 
forms of the V(TCNE)2 compound as the highly magnetic (HM) and the low 
magnetic(LM) ones. 

For more than one decade the amorphousity of the compound of inter- 
est did not allow anyone to make whatever definitive conclusion concerning 
its structure. The critical breakthrough became possible with the recent 
work [8] where the authors were able to establish the structure of the Fe^"*" 
analog (presented in Fig. 1) of the V(TCNE)2 compound using Rietveld re- 
finement of the synchrotron powder diffraction data and to reveal its most 
remarkable features: the presence of the dimer form of the TCNE~ radical- 
anion: [TCNE]2~ = C4(CN)|~ playing an important role in shaping the loose 
three-dimensional structure and assuring as well the three-dimensional char- 
acter of magnetic interactions in the system. 

In our previous paper Ref. [10] we were able to demonstrate that the 
experimental structure represented in Fig. 1 can be easily related with one 
represented in Fig. 2a proposed yet in Ref. [9] in order to conform with the 
magnetic data on the V-TCNE compound available that time. That latter 
represented a simple cubic lattice with the vertices occupied by the vana- 
dium ions. Two a6-faces of each cube remain empty whereas four others are 
filled by the TCNE units forming channels extended in the c-direction. In 
the structure presented in Fig. 2a the C=C bonds of the TCNE units are 
located in the same plane (as well perpendicular to the c-axis of the lattice) 
so that they are orthogonal to each other. It is not the unique possibility. 
An alternative structure differing from that of Fig. 2a by rotating the TCNE 
units placed in the 6c-faces by 90° in their respective planes is also possible. 
The result of such a rotation is presented in Fig. 2b. In this structure the 
central C=C bonds of the TCNE units previously orthogonal, but 

now they lay in orthogonal planes so that the axes of these bonds do not 
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intersect rather cross each other. This structure has been called the "prin- 
cipal" structure in Ref. [10]. The principal structure Fig. 2b can be easily 
put in the relation with the experimental one. If in the quadrupled unit 
cell 2a, 2b, c of the principal structure four TCNE units extended in the b- 
dircction are allowed to pairwisely rotate towards each other so that single 
C-C bonds can form between respective cthylcnic carbon atoms thus yielding 
the [TCNE]2"=C4(CN)8" dimers and the V-TCNE sheets originally laying 
in the ac-planes are accordingly ruffled one finally arrives to the experimental 
structure of Fe(TCNE)2 (Fig. 1). Intermediate structures along this hypo- 
thetical reaction path are presented in Fig. 3. For this sequence of structures 
we performed in Ref. [10] the LSDA + U calculations with use of the VASP 
program suite Ref. [11] of the respective electronic structures and energies. 
Specifically, the PAW potential has been used and the values of the U param- 
eter for the V, N, and C atoms were respectively taken as 5.344, 4.840, and 
4.428 eV. The calculations have been performed as follows: first the principal 
structure with the equalized lattice parameters have been quadrupled. This 
structure has been connected with the experimental structure by a straight 
line in the configuration space. Further details of the numerical treatment 
can be found in Ref. [10]. Its results are reproduced in Fig. 4. 

The overall densities of states in two spin channels as given in Fig. 4 are 
obviously difficult to understand. One can only see some spin polarization 
of the upper filled bands as well as an evolution of a noticeable density of 
states near the Fermi level present in both spin channels in the initial state to 
the final state where the expected spin-polarized structure can be recognized 
which docs not show any DoS in either of the spin channels at the Fermi level. 
In this paper we present a detailed analysis of the results of our numerical 
studies Ref. [10] of the models of V-TCNE room temperature organometallic 
ferrimagnet and expressed them in terms of the effective spin-Hamiltonian 
for a selection of interacting atomic/molecular states. The proposed model is 
then applied to analysis of a wider collection of experimental data available 
for this fascinating object and its Fe analog. 

2 Detailed analysis of DoS 

In order to analyse the details of the evolution of the DoS features along the 
path going from quadrupled principal structure to the experimental one we 
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performed a detailed study of projections of the DoS. In the initial structure 
the three bands in the —9^ — 11 eV range are predominantly C-bands. Other 
bands are the well hybridized C-N-bands except the "spin-left" band right 
below the Fermi level, which is predominantly contributed by the vanadium 
states. As one can see from Fig. 6 these are the rf-states of vanadium atoms. 

Some vanadium stemming DoS in the range —8 10 eV are the s- and 

p-states of vanadium involved in the bonding with the nitrogen atoms. This 
can be seen on the corresponding N-projection of DoS in Fig. 7 from which 
one can deduce that the s- and p-states of vanadium hybridize predominantly 
with the N-states coming from the V-TCNE layers (see below). In order to 
deepen our understanding of the DoS presented in Fig. 5 we notice that 
even in that complex system the atoms involved are relatively easily classi- 
fied into types. One can distinguish transition metal ions whose d-density is 
expected to contribute significantly to the spin polarized bands, the nitrogen 
atoms where one can expect significant DoS changes due to expected break of 
the pairs of excessive V-N bonds while going from the quadrupled principal 
structure to the experimental one. Analogously along the same route one 
can expect remarkable variations in the projection of the DoS to the carbon 
atoms forming the C-C bonds in the [TCNE]|~ units and to the carbon and 
nitrogen atoms in the ruffled V-TCNE layers. Other features of the DoS, 
however, are not expected to significantly vary along the path. Considering 
the evolution of the above-mentioned projections of the DoS along the path 
Fig. 3 as represented in Figs. 4-6 one can see the expected reconstruc- 
tion of the projected DoS's. For example, almost nothing happens to the 
rf-bands along this path. They remain rigid triply degenerate ones and do 
not change their position relative to the Fermi level despite the fact that the 
coordination number of vanadium ions changes from eight in the quadrupled 
principal structure to six in the experimental one. This agrees with the nu- 
merical result concerning the distribution of spin polarization in the direct 
space: for all structures depicted in Fig. 3 the magnetic moments residing 
in the rf-shells of vanadium ions range from 2.615 for the first (quadrupled 
principal) structure to 2.582 for the last (experimental) one i.e. are almost 
constant. The significant variation of the total magnetic moment along the 
"reaction path" observed in our numerical experiment is to be almost com- 
pletely attributed to that of the magnetic moments residing in the "organic" 
part of the organometallic magnet. This is precisely what one should expect 
within the general picture including the formation of the [TCNE]2~ dimers. 
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On the other hand the stabihty of the rf-bands along the path can be only 
understood if one assumes that the break of two V-N bonds at each atom 
is at least partially compensated by shortening (strengthening) of other two 
bonds extended in the 6-direction. 

By contrast the projections of DoS to the characteristic organogenic atoms 
significantly modify along the "reaction path". For example, the dominant 
part of the N-DoS is concentrated in two broad bands at ca. —5 and ca. 
—6-. — 6.5 eV. The lower of the two is almost equally contributed by the 
nitrogen atoms of all three types and its position is stable throughout the 
path. The same applies to somewhat smaller contribution from the bonding 
N-atoms extended in the 6-direction to this band. Their contribution to the 
upper of the two mentioned bands is noticeably spin-polarized. The density 
from the N atoms bearing the dangling lone pairs fairly manifests itself in 
the same band. Weak features in the N-DoS can be observed in the 4 eV 
wide range right below the Fermi level. The peak right below the Fermi level 
in the "right-spin" channel in the initial structure is equally contributed by 
the "bonding", "to be dangling", and "magnetic" nitrogens. (In general the 
DoS projections for the "bonding" and "to be dangling" nitrogens coincide 
in the initial structure). In the final structure this peak splits so that its 
upper part ("right-spin" channel right below the Fermi level) is contributed 
by the "magnetic" nitrogens, whereas two peaks next to it in the bottom 
direction is contributed by the "bonding" nitrogens. The contribution from 
the "dangling" nitrogens contributes to the upper part of the wide band at 
ca. —3.5 -= — 5.5 eV. 

Of particular interest is the evolution of the C-projection of the DoS. That 
on the ethylenic carbon atoms is most interesting. The cthylcnic carbons 
can be further subdivided in two categories: those in the planes extended 
in the a and c horizontal directions and those extended in the vertical (6) 
direction and further involved in the rotation yielding the [TCNEjg" units. 
At the initial stage both projected DoS are significantly polarized and both 
noticeably contribute to the "right-spin" density right below the Fermi level 
(the DoS projection to the to be bonding atoms is not seen in the leftmost 
graph of Fig. 6 since at this point they are degenerate with the "magnetic" 
projection and are masked by these latter). Prom further graphs of the DoS 
projected to the atoms forming the emerging C-C bonds one can conclude 
that such bonds represented in our study by two spin sub-bands for the 
left and right spin channels completely develop at a pretty late stage of 
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the hypothetical transition: in the middle of the path the corresponding 
peaks in the density projections can be yet clearly seen in both spin channels 
near to the spin-polarized DoS of the "magnetic" atoms, although the peaks 
corresponding to the emerging bonds are not polarized. By contrast the 
projection of DoS on the magnetic C-atoms in the horizontal planes develop 
two sub bands right above and below the Fermi level of which the lower 
one (right-spin) is completely occupied by electrons with the spin projection 
opposite to that of the electrons occupying the d-subband. This is in the 
fair agreement with the assumption concerning the nature of the subbands 
located in the vicinity of the Fermi level made in our previous paper Ref. [10]. 

3 Model Hamiltonians for V-TCNE system 
3.1 Model Orbital Hamiltonian 

Fascinating properties of the V(TCNE)2 magnet call not only for numerical 
modelling, but also for some qualitative picture. However, for the final ("ex- 
perimental") structure all the DoS projections manifest themselves as very 
narrow bands. This indirectly indicates that the band picture used through- 
out the calculations is not completely adequate and that an adequate model 
must be given in terms of an effective Hamiltonian representing the electronic 
structure of the HM V(TCNE)2 magnet (in its "experimental" structure) in 
terms of some objects local in the direct space e.g. local spins similar to that 
proposed yet in Ref. [9]. As in the case of band models the most important 
one-electron states to be included are these contributing to the energy bands 
in the vicinity of the Fermi level. Based on our analysis of projected DoS 
performed in the previous Section we can conclude that the states of the 
[TCNE]|~ units contribute to the bands far away from the Fermi level. The 
DoS related to this unit goes away from the Fermi level along the "reaction 
path". Thus the observed electronic structure is primary one of the indi- 
vidual (ruffled) V-TCNE layer extended in the ac-plane. For constructing 
the Hamiltonian for this layer one can employ the unit cell of the principal 
model dropping from it the TCNE unit extended in the 6-direction (and fi- 
nally engaged in formation of the [TCNE]2~= C4(CN)|~ dimers). According 
to Ref. [9] on the vanadium sites it suffice to consider only the d-shells of the 
metal ions which is also confirmed by our current analysis. The overlap of the 
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(i-shell of the metal ions with the cr-orbitals of the TCNE's (inchiding those 
imphed in the model) ensures the standard two-over-three splitting of the d- 
shell characteristic for the octahedral environment. In the case of vanadium 
three unpaired electrons in the rf-shell occupy respectively three orbitals in 
the t2g-va.amlo\(\. The basis orbitals d^-y^ d^z, and dy^ can be characterized by 
the normal to the plane in which each of the orbitals lays - (, r], and ^ - are 
subsequently used in the notation. For the donor sites in the ruffled planes 
the bsg (it*) LUMOs of the TCNE (singly occupied in the radical anion) are 
included. In each such a layer each metal ion is surrounded (coordinated) 
by four TCNE units which are in their turn coordinated to (surrounded by) 
four metal atoms. In this case the layer unit cell composition V:TCNE is 
1:1. 

The model Hamiltonian for the V(TCNE)2 magnet, formulating the above 
ideas, has the general form: 

H = J2 (^d(r) + Hair) + H^air) + H^ir)) (1) 

r 

The contributions to it are the following. Operator Ha{r) describes electrons 
in the acceptor orbital of the TCNE" radical-anion in the r-th unit cell: 

Hair) = 

Symbol a^^{ara) is the operator creating (annihilating) an electron with the 
spin projection a on the acceptor orbital of the TCNE molecules in the r-th 
unit cell. In eq. (2) the first term is the energy of attraction of an electron to 
the core of TCNE - the orbital energy of the 63^ (tt*) LUMO shifted by the 
electrostatic field induced by the entire crystal environment. The second term 
in eq. (2) is the Hubbard one, effectively describing the Coulomb repulsion 
of electrons with opposite spin projections eventually occupying the same 
acceptor orbital. 

The operator H^ir) describes electrons in the i2g-subshell of the d-shell 



L-ia "'ara 
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of the vanadium ion in the r-th unit ceU: 



+ {Udd + '^Jdd){nCv\.n(^Tt + flr^Tiflr^r^ + n^r4.%t)] + 
, {Udd+ Jdd/'^) sr^i ^ . , . . , - - \ , /Q\ 



n. 




a 



In eq. (3) n^ro- are the operators of the number of electrons with the spin 
projection a on the dxy, dy^, and d^^ orbitals of the vanadium ion in the r-th 
unit celL The spin operators and spin-operator product terms are defined by 
the well-known relations: 



where the symbols ilra) represent the operators creating (annihilating) 
an electron with the spin projection a on the dxy, dyz, and dxz orbitals of the 
vanadium ion in the r-th unit cell. 

The first row in the above operator describes the attraction of electrons in 
the d-orbitals to the cores of vanadium ions (shifted by the electrostatic field 
of the rest of the crystal). Two further rows describe the spin-symmetric 
part of the Coulomb interaction of electrons in the (i-shell. The last row 
describes the spin dependent part of the Coulomb interaction of electrons 
in the (i-shell (exchange). It is ultimately responsible for the Hund's rule in 
atoms and for the high spin of the ground state of electrons in the d-shell. 
The contributions to the Hamiltonian eq. (1) described so far model isolated 
local states important for the crystal description. The magnetic order can 
only be possible due to various interaction terms. Operator i/(ia(r) describes 
the electron hopping between the (i-states of vanadium ions and the acceptor 
states. The dxy-state represented by the (^^ (Cra) operators being of the 
(approximate) cr-symmetry with respect to the ac plane (the ruffling of the V- 
TCNE plane is neglected) has no overlap with the LUMO's of TCNE's which 
are (again approximately) of the vr-symmetry with respect to the same plane. 
Two others {d^z- and d^^-states represented respectively by the {Vra) and 



l/2(iS'^j.iS'_,/j. + iS'„,/i.»S'-,j.) + S^j.S, 
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^rcr (Crcr) operatois) overlap with the LUMOs of two (different) neighbor 
TCNE units each. The phase relations between the orbitals involved in 
the model lead to such a distribution of signs at the one-electron hopping 
parameters that the hopping operator acquires the form: 



where the parameter t^a > describes the magnitude of the hopping between 
the acceptor state and the neighbor (i-state. 

The sum of the above contributions to the effective Hamiltonian in fact 
form that for an isolated V-TCNE layer. In the "experimental" structure 
the diamagnetic C4(CN)8~ units seem to effectively isolate the V(TCNE) 
sheets from each other. Nevertheless, one should assume that certain indirect 
interaction between the d-states in the 6-direction is possible through the 
mediation of the [TCNE]2~ units. It was proposed in Ref. [10] to use an 
effective hopping similar to eq. (4) . Since it is any way an effective interaction 
it can be chosen in a way which fits better to the method the system is treated. 
For this reason we postpone the discussion of this term. 

In our previous paper Ref. [10] we considered the band model of the V- 
TCNE organometallic magnet as derived from the orbital Hamiltonian eqs. 
(1) - (4) and employed them for analysis of results of our numerical experi- 
ments performed with use of the VASP package. These latter are, however, in 
a kind of fundamental contradiction with the physics of the system at hand. 
This manifests itself in the very narrow bands coming out of calculation, as 
we already mentioned. The reason is that the hopping parameter tda entering 
eq. (4) which are generally responsible for extension of one-electron states 
over the crystal (band formation) and which are proportional to the overlap 
between the orbitals represent the smallest energy scale in the system. Gen- 
erally it leads to a break of the delocalized (band) picture and makes a local 
description to be more adequate. The latter can be sequentially derived by 
treating perturbatively the hopping operator eq. (4). It yields the effective 
Hamiltonian of the Heisenberg form in terms of the spins of electrons occu- 
pying the local states (orbitals) involved. Its parameters are estimated in 
Appendix A. 



Hdair) = -tda l^ra («ra + 



O-r+a+ca) ~ Vta ('^r+atr + O'r+ca) + h.C. 



(4) 
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The overall result comes out as a spin Hamiltonian of the form: 



r 

+2i^^do S S(^r [par + + 



(5) 



which describes effective magnetic interactions in an isolated layer. It must 
be complemented by interlayer interactions. If the vanadium ions in adjacent 
layers are coupled by an an effective hopping an analogous perturbative pro- 
cedure results in the antiferromagnetic sign of the effective magnetic interac- 
tion. This contradicts to the existence of the nonzero overall magnetization in 
the V-TCNE magnets below the critical temperature (with the antiferromag- 
netic interlayer coupling the magnetization of one layer would be cancelled 
by that of another). For that reason we have to supplement the Hamiltonian 
eq. (5) by an effective interlayer interaction with the ferromagnetic sign of 
the corresponding exchange parameter. It cannot directly come from any 
perturbative treatment of the hopping. By contrast some mechanism of fer- 
romagnetic coupling described e.g. in Refs. [12] or [13] and implemented in 
papers [14, 15] devoted to the exchange in metallocene based organometallic 
magnets (Miller-Epstein magnets) acting through the [TCNE]2~ units might 
be expected. Indeed as one can see from the Fig. 7 despite the fact the 
the states located in the [TCNEjg" units are pulled up and down from the 
vicinity of the Fermi level, some spin polarization of these bands particu- 
larly of those which are contributed by the "bonding" nitrogens indicates 
the involvement of the [TCNE]2^ units in transfer of magnetic interactions 
between the d-shells in the 6-direction (between the layers). 

With this caveat the spin Hamiltonian written in terms of "true" elec- 
tronic spins is the following: 

Hspin — ~4:Jfj^fj^J2{Si^rS^r~^S;^j.S^j.-\-S^j^Srir)'\' 

r 

+ '^Kda Til S^r (^Sar + >S'ar+a+c) + 'S'^r (^'S'ar+a + 'S'ar+c) (6) 
+ 2/^^^ ('S'^r + Sjjr + Si^r^ ('S'^r+b + 'S'jjr+b + 'S'^r+b) 

The Hamiltonian eq. (6) is not a standard Heisenberg Hamiltonian usually 
used to describe magnetic properties of insulators. This latter is written in 
terms of the effective local spins residing at each atomic magnetic center. 
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In our case the vanadium ions represetn such nontrivial magnetic centers 
bearing effective spins in the ci-shells: 



SdT — '^Sjr- (7) 

According to Ref. [16] using the effective spins is, however, an approximation 
since the transition from the representation of the effective Hamiltonian in 
terms of the of individual electronic spins-| eq. (6) which can be sequentially 
derived from the model orbital Hamiltonian eqs. (1) - (4) by perturbative 
treatment of the hopping term eq. (4) to the phenomenological Hamiltonian 
eq. (8) operating with the effective spins eq. (7) is only possible if the 
exchange interactions of all individual spins in one magnetic center (in our 
case - the V ion) with those in the other magnetic center (in our case the 
effective spin in TCNE~ coincides with the individual one) are equal. This is 
obviously not the case since the electronic spin S^r in the d-shell to the first 
approximation does not interact with the spin residing in any of acceptor 
orbitals. 

This generally poses the problem since in the Hamiltonian eq. (6) at 
least the exchange parameters Jdd and Kda can be independently determined 
respectively by eq. (18) and atomic spectra, whereas the parameter J± in 
eq. (8) remains completely empirical quantity. The fact that Jdd 3> Kda in 
eq. (6) allows to approximately replace the spins of separate electrons in the 
r-th d-shell by the operator of the total spin of the respective (i-shell. Further 
details of this transition are given in Appendix B. 

The phenomenological Hamiltonian written in terms of the effective spins 
eq. (7) to be used for modeling the entire crystal is: 



Hphen — J\\Y1 Sdr Sar + 'S'ar+a + 'S'ar+c + Sar+a+c 
r L 

+ J_L SdrSdr+h 



+ 

(8) 



In the next Section we apply it to analysis of magnetic properties of the HM 
V(TCNE)2 material. 
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4 Magnetic properties of V(TCNE)2 as inter- 
preted with use of phenomenological Hamil- 
tonian 

The magnetic and thermodynamic properties of the HM V(TCNE)2 material 
must be derived from the phenomenological Hamiltonian eq. (8). Two types 
of data will be of interest for us: the critical temperature of transition into 
magnetically ordered state and the temperature dependence of spontaneous 
magnetization. The mean field estimates of the critical temperature used so 
far in the literature lack the account of structural information. It is impor- 
tant to realize that the quantities of interest are sensitive to these details as 
represented in the respective Hamiltonians. 

The commonly used (see e.g. Ref. [6]) symmetric mean field formula: 

er = '-^Sa {Sa + 1) 

(where T = 9/kB) ignores the acceptor spins and sets z to be the number 
of indirectly neighboring magnetic metal ions. It yields the quadratic de- 
pendence of the critical temperature on the spin of the metal ion. On the 
other hand according to Ref. [8] the critical temperature for the magneti- 
cally ordered state of a material with two types of spins {Sd and Sa{= 1/2)) 
is described by the mean field formula: 



9 



MF 
N 



tMF 



JSa{Sa+l)Sa{Sa + l) 



(The factor of two is dropped here to get the formula to conform with the 
Hamiltonian definition accepted in the present paper). This formula appears 
as a zero interlayer coupling limit of the mean field expression for the Neel 
temperature in a ferrimagnet eq. (35) derived in Appendix C. In Ref. [8] 
it had been applied to the Fe(TCNE)2 compound for which the structure 
measurements have been performed there. It, however, brings up two com- 
plications - one theoretical and another experimental. From the experimental 
point of view we notice that the above formula as well as formula eq. (35) is 
effectively linear in Sd rather than quadratic (in the high anisotropy limit). 
For that reason even the mean field estimates of the exchange parameters 
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as given in Table 1 must be reconsidered since the latter had been obtained 
with use of the quadratic dependence. As one can see from the structure the 
choice of Zad — Z^a — 4 yields the estimate of the mean field exchange param- 
eter for the Fe(TCNE)2 compound of J^^{Fe) = 43 K and for the V(TCNE)2 
compound - J^^(V) = 183 K. Both values are significantly larger than those 
given in Table 1, but it is remarkable that the difference between them (to 
be explained) reduces from the factor of larger than five to that of 4.25. We 
notice that following the assumption of Ref. [24] and using Zm — 6 for the 
V(TCNE)2 compound further reduces the difference and the value of the 
exchange parameter for the latter, but from our point of view it cannot be 
substantiated within the scope of the model considered in the present paper. 

From the theoretical point of view, even the improved molecular field 
expression eq. (35) has that disadvantage that it predicts a nonvanishing 
ordering temperature for J± = 0. It is obviously wrong and such an es- 
timate is not acceptable in the context where a strong anisotropy might be 
expected on the structure basis. As it has been mentioned the model must be 
complemented by the interlayer interactions between the effective spins 3/2 
on the vanadium sites mediated by the diamagnetic (closed shell) [TCNE]^~ 
units. Remarkably enough the sign of this interaction must be ferromag- 
netic (the magnetic moments residing in the layers must be pointing in the 
same direction) to ensure the existence of the net spontaneous magnetiza- 
tion in the three-dimensional sample, although in general one has to expect 
antiferromagnetic sign of such an interaction [12] (see below). 

In a presumably rather anisotropic situation brought by tentative differ- 
ence in mechanisms of the intralayer and interlayer interactions (respectively 
"antiferromagnetic kinetic exchange" for the intralayer interaction and the 
"ferromagnetic superexchange" for the interlayer one) the critical tempera- 
ture has to be estimated from the spin-wave treatment taking an adequate 
care about the anisotropy of the effective spin-spin interaction and at least 
providing a correct asymptotic value of the critical temperature for the van- 
ishing interlayer coupling J±. This is done by the formula 



expressing the Bloch Ta law for the temperature dependence of the spon- 
taneous magnetization as derived in Appendix D with the critical (Neel) 
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temperature given by: 



= i!,,^ ^ ^ ^{Sd + S^fSjSl [Jf |JJ]^ (10) 

- the 3D structure-specific relation of the effective exchange interaction to 
the Neel temperature. 

The assumptions used in Appendix D for deriving the formula eq. (10) are 
not satisfied in strongly anisotropic systems where | J_l| <C Jy. In this limit 
the Neel temperature is to be determind from the transcendental equation: 



1 Sd — Sn ^> 



47r SdSa J\\ Sd + Sa 

and the magnetization is given by: 




1 Sd — Sa 1 , f Sd — Sa 

log 7^ 



AtT SdSa Sd + SaJ\\ 




(11) 



(12) 



provided 9 is close enough to ^n- 

Neither three-dimensional (3D) or two-dimensional (2D) estimates of the 
Neel temperature eqs. (10) and (11), respectively, permits to determine 
the longitudinal and transversal interactions independently and to estab- 
lish by t his th e amount of anisotropy. The effective exchange interaction 
JeS^ = \fjfji. ^ derived from eq. (10) and the experimental Neel tempera- 
ture for the HM V(TCNE)2 compound amounts Jf^(V) = 36 K. This is due 
rather large numerical value of the transition coefficient in eq. (10) coupling 
the effective exchange interaction with the Nccl temperature (11.376 for the 
structure depicted in Fig. 1 and Sd = \; Sa = \)- This result is in a general 
agreement with the result of Ref. [18] which yields the corresponding coeffi- 
cient to be 9.937 for a simple cubic ferrimagnet with the same values of the 
effective spins. (It is not clear how the estimate of ca. 100 K for Jcf^(V) 
is obtained in Ref. [23] since is also based on assumption of a simple cubic 
lattice magnetic structure, but apparently uses some different coefficient). 
It also stresses the different character of averaging of intralayer and inter- 
layer exchange parameters in the mean field (arithmetic mean) and in the 
spin-wave (geometric mean) approximations. At the high anisotropics the ge- 
ometric mean provides much stronger dependence of the effective exchange 
on the interlayer exchange than the arithmetic mean. 
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Whatever value of leaves a wide range of possibilities since each pair 
of values of J± and Jy yielding the above value of ^eff^(V) conforms with the 
experimental data on magnetization. It has to be realized, however, that us- 

3 

ing the Bloch T2 law for the magnetization in the entire temperature range 
below Tjv is an extrapolation of the data obtained at low temperatures. In 
order to check its validity in a wider temperature range we notice that ac- 
cording to it the magnetization depletion at the an intermediate temperature 
(225 K) amounts the factor of 0.592 which looks out to be in an acceptable 
agreement with experiment which shows the magnetization depletion by a 
factor of ca. 0.6 at this temperature as compared to that at T = 0. 

In the HM V(TCNE)2 case the measured magnetization values are avail- 
able up to 300 K. The Bloch Tz-law for the temperature dependence of 
spontaneous magnetization results in a simple formula for the slope of the 
magnetization vs. temperature in the Neel point: 

-l^ (13) 

As one can derive from the magnetization data on the HM compound given 
in Ref. [6] the slope of the magnetization of the HM material in the Neel point 

amounts —1.4—- which is a fair extrapolation of the last measured points. 

It suggests the 3D regime for the HM material in the entire temperature 
interval up to the Neel point. Nevertheless the possibility of transition to 
the 2D regime at T > 300 K cannot be a priori excluded. Then for the 2D 
regime the slope of magnetization vs. temperature in the Neel point is: 



ks Sd — Sa 1 / T , 1 f Sd — Sa ksTi 



When combined with eq. (H) it yields: 

-Mo ( ^d-Sg 1 _^ ^\ ^^^^ 
\47rJ|| SdSa Sd + Sa T^j 

Employing the value of the slope extracted from Fig. 1 of Ref. [6] we de- 
rive /cbTn — 2.47rJ||, and inserting experimental (extrapolated) value of Tn 
yields immediately .J\\ = 54 K and together with the value of .J^^ extracted 
from low temperature data allows to estimate anisotropy to be {J\\/ J±_ ~ 3). 
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At the above intermediate temperature (225 K) the 2D estimate with this 
anisotropy shows the depletion of magnetization to be 0.595 of the maximal 
value at K as well in a perfect agreement with experiment, which shows 
that the available data on the temperature dependence of magnetization in 
HM V(TCNE)2 do not allow to distinguish between the 3D and 2D regimes. 
It must be admitted that in general the above vahic of anisotropy is not large 
enough (~ 3) for the 2D regime to install. This analysis, however, allows us 
to set bounds for the value J\\ in the V(TCNE)2 compound as derived from 
the spin- wave treatment: it appears that this compound resides in the 3D 
regime so that the above value of anisotropy must be considered as a maximal 
possible in this material. Otherwise even higher Neel temperatures (although 
not accessible expreimentally due to material's decomposition) still conform- 
ing to the applicability conditions {'2J±Sj < 6* < AJ^^SaSa) Ref. [25] of the 
logarithmic formula eq. (11) should have to be admitted. 

Applying analogous treatment to the Fe(TCNE)2 compound for which 
the structure presented in Fig. 1 is experimentally established yields the 
following: J^g^ = 9.5 K (the coefficient of 12.914 coming from eq. (10) with 
Sd — "2 is used). With the anisotropy of 2.5^ = 15.625 the 2D estimate of the 
Neel temperature is 124 K again in a fair agreement with the experiment. 
The value of Jy is then 24 K. This turns out to be not that much different 
from the upper boundary for the same quantity for the V(TCNE)2 compound 
yielding the ratio of the intralayer exchange parameters for the two materials 
of maximum only two, instead of 4 5 stipulated by the mean field estimates, 
and only 1.5 if the isotropic regime is accepted for the V(TCNE)2 compound. 

This latter value can be fairly explained by addressing the formulae of 
Appendix A and the spectroscopic data. From eq. (18) it follows that the 
ratio J|| (V) / Jy (Fe) of the intralayer parameters for the vanadium and iron 
compounds is that of the squared hopping parameters tda- (We assume here 
that due to similarity of the environment in these two compounds the en- 
ergy denominators in eq. (18) given by eq. (19) are the same for the both 
compounds since the values of ionisation potentials of the V^^ and Fe^+ ions 
which are respectively 29.55 and 30.90 eV as coming from Ref. [26] and sim- 
ilarly close estimates for the electron affinities for these ions) . According to 
suggestion by [27] thoroughly tested numerically in Refs. [28-30] the amounts 
of the crystal field splitting in the coordination compounds are proportional 
to analogous expressions: squares of the hopping parameters divided by some 
(other) energy denominators, which are, however, also approximately equal 
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in similar compounds. Thus the proportion holds: 



J||(V) _ lODg(V) 
J||(Fe) ~ 10L>g(Fe) 

for pairs of similar complexes in each side of the proportion. For the right side 
of the proportion we find with the values of Refs. [31,32] 10Dq{V) = 14700 
cm~^, lODg(Fe) = 10900 cm~^ to be 1.35 for the hexacoordinate octahedral 
complexes with acetonitile, which basically explains the above ratio 1.5 of 
the intralayer exchange parameters. Of course, the significant difference of 
anisotropies in these two materials remains to be understood. 

This can be tentatively done with use of the Goodenough-Kanamori 
rules [12]. Indeed, the difference in the interlayer interactions requiring an 
explanation is too large, so that probably a qualitative distinction between 
the two materials is responsible for it. As we mentioned above the simplis- 
tic application of the Goodenough-Kanamori rules in the present situation 
yields an antiferromagnetic sign of the interlayer interaction (the situation 
falls into the Goodenough-Kanamori cation-anion-cation category in the 180° 
geometry). Thus the observed ferromagnetic sign of the interlayer interac- 
tion appears as a result of ferromagnetic contributions of the higher order. 
Such contributions depend qualitatively on the possibility to take advantage 
of the intrashell ferromagnetic interactions which in their turn depend on the 
occupancies of the atomic orbitals in the d-shells of the interacting transition 
metal cations. (Importance of such terms in the context of organometallic 
magnets had been stressed in Refs. [14,15]). It can be easily understood that 
the conditions for appearance of the compensating ferromagnetic terms are 
very much different for the and Fe^+ ions. Indeed, in the case of the 
V^"^ ion two d-orbitals remain empty and can participate in the one-electron 
transfer coupled with the intrashell exchange eq.(21) compensating other- 
wise dominating antiferromagnetic kinetic exchange. In the case of the Fe^"*" 
ion only one doubly occupied d-orbital can take part in a similar process, so 
that one can expect that the compensating contribution will be significantly 
weaker in the case of the Fe(TCNE)2 compound eventually leading to much 
weaker overall ferromagnetic interlayer interaction, than in the case of the 
V(TCNE)2 compound. 
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5 Discussion 



In the present Section we apply the models proposed above to analysis of 
expcrimantal data available for the HM V(TCNE)2 and Fe(TCNE)2 com- 
pounds. 

First of all we notice that the spin polarization per unit cell (number 
of electrons with spin up minus that with spin down) which can be related 
with observed magnetization per formula unit. We see that the calculation 
performed for V(TCNE)2 at the experimental structure of Fe(TCNE)2 de- 
picted on Fig. 1 shows the spin polarization of ca. 8 spins- 1/2 per unit cell 
corresponding to two netto unpaired electrons per formula unit which is in a 
fair agreement with the magnetization measured in the HM V(TCNE)2 com- 
pound. On the other hand the LM V(TCNE)2 material manifests a weaker 
saturation magnetization, namely corresponding to ca. one netto unpaired 
electron per formula unit. This allows to think about certain differences in 
the structures of two materials. Nevertheless, both experimentally observed 
values (ca. 10-10^ emu-Oe-mol~^ and 6-10^ emu-Oe-mol~^, respectively) both 
deviate from the theoretical values of 11.2 and 5.6 giving the magnetization 
produced by the integer number of netto spin-polarized electrons in an as- 
sumption of the Landc factor being equal to 2. 

When trying to extend the model Ref. [8] of the Fe(TCNE)2 compound 
to analysis of the V (TONE) 2 compound the authors Ref. [8] argued that 
the vanadium compound must have some structure different from the iron 
one since the saturation magnetization in it is lower and approximately cor- 
responds to two spins 1/2 compensating (interacting antiferromagnetically 
with) one spin 3/2 per formula unit. From this observation the authors 
of Ref. [8] conclude that the interlayer interactions must be mediated by 
//4-TCNE radical-anions, as it has been suggested yet in [9], rather by the 
[TCNE]2^ dimers. This argument applies of course only to the LM form 
of the V(TCNE)2 material since for its HM form our numerical experiment 
shows that for the experimental structure of the Fe(TCNE)2 compound the 
calculated magnetization fairly corresponds to the experimental value ob- 
tained on the HM V (TONE) 2 material. Incidentally, the magnetization val- 
ues obtained numerically at intermediate structures on the "reaction path" 
depicted on Fig. 3 allows us to assume that some similar structures obtained 
from structures of Fig. 2 by rotations of some TONE units may present 
in the LM V-TCNE material. This view had found a recent support from 
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the computational side in Ref. [33] where a structure for the LM form of 
V(TCNE)2 has been proposed. It would be fair to say (although it is not 
said in Ref. [33]) that this structure as well descends from the structures of 
Refs. [9, 10]. Specifically, in order to obtain the structure of Ref. [33] one has 
to rotate the TCNE molecule laying in the 6c-face of either of these structures 
depicted in Fig. 2 in each of the unit cells around the diagonal of the face (or 
around an axis going through the pair of trans-nitrogen atoms of that TCNE 
unit) by ca. 90° so that two other N-atoms of each rotating TCNE unit go 
out of coordination with the V ions. Such a structure corresponds as that of 
Ref. [9] and the principal one (see above) of Ref. [10] to two TCNE^ units 
per unit cell each bearing one unpaired electron and thus expectedly yield 
the overall magnetization corresponding to one unpaired electron per formula 
unit. For such a structure the magnetic interaction parameters obtained in 
Ref. [33] arc almost isotropic (J_l = 720 K and = 690 K) which is also 
not surprising since the character of interactions between the V rf-shells and 
LUMO's of the TCNE^ units arc fairly the same in either direction. The 
numerical values of the exchange parameters obtained in Ref. [33] are for 
sure considerable overestimates of the true ones since the Neel temperature 
derived from them either by the mean field or spin-wave methods exceeds the 
experimental value by orders of magnitude. Nevertheless, applying the spin- 
wave theory similar to that described in Section D results in the estimate 
for /cbTn = 9.14Jgg^ which yields the numerical value of for the LM 
form of V(TCNE)2 material of 45 K in fair agreement with the similar above 
estimates for the intralayer effective exchange parameters. It must also be 
admitted that the spin-wave treatment of the model of Ref. [33] leaves the 
question of the reason of complete disagreement of the temperature depen- 
dence of the magnetization in the LM compound as given in Ref. [6] with 
the Bloch law which should be expected for almost isotropic ferrimagnet 
unanswered. 



6 Conclusion 

In the present paper we performed detailed analysis of our numerical results 
concerning thinkable structure of room-temperature organometallic magnet 
V(TCNE)2 as manifested in the corresponding projections of DoS. Similar 
analysis of projected DoS for a sequence of structures leading to the ten- 
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tative experimental structure of V(TCNE)2 is performed as well. Model 
spin Hamiltonian is developed for analysis and interpretation of numerical 
results and experimental data. Analysis of magnetic data in terms of the 
approximate models derived from the phenomenological Hamiltonian is per- 
formed. A remarkable correspondence between experimental (structural and 
magnetic) data on V(TCNE).^. • y solvent and numerical model has been ob- 
served previously: magnetization corresponding to two unpaired electrons 
per formula unit in fair agreement with experiment on HM V-TCNE ma- 
terial derived from V(C0)6 by CVD technique is obtained numerically for 
V(TCNE)2 taken in the relaxed experimental Fe(TCNE)2 geometry Ref. [10]. 
Now it is complemented by the detailed analysis of the magnon spectrum of 
this model. The possible transition between the low-temperature 3D and the 
high-temperature 2D regimes is discussed. Estimates of parameters of the 
proposed spin- Hamiltonian as treated in the spin-wave approximation are 
derived from the experimental data on the Neel temperature and the tem- 
perature dependence of magnetization. The differences in magnetic behavior 
of probably isostructural HM V(TCNE)2 and Fe(TCNE)2 are tentatively ex- 
plained. 
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A Spin Hamiltonian as derived from orbital 
model Hamiltonian 



Simple estimates of the parameters of the spin Hamiltonian eq, (6) can be 
based on considerations dating back to Ref. [12] as specified for the current 
situation. For the idealized geometry of M(TCNE)2 compounds where the M- 
TCNE layers are assumed to be planar the one-electron hopping parameters 
tda are those between the b^g (tt*) singly occupied orbital of the TCNE~ 
radial-anion and one of the d-orbitals of the 7r-symmetry with respect to the 
abovementioned plane (d^z and dyz). For a pair of the TCNE~ and V^"*" ions 
described by the Hamiltonians, eqs. (2) and (3), respectively, the energy of 
the bare ground state reads: 

Eo = -aa - Sad + SUdd - 2l/AJdd (16) 

and does not depend on the way the spins in these sites are coupled. The one- 
electron hopping couples them with two states with one electron transferred 

between the two sites (from one of the rf-states to the a-state and from the 
a-statc to one of the d-states). The energies of the charge transfer states are 
respectively: 

Ed^a = -"i-OLa - 2«d + Uaa + U dd " 3 Jdd/2 

E^^d = -Aad + 6Udd + 29/AJd4 

They both correspond to the lower spin of the state of the two ions which 
results, as usual, to a Heisenberg-type interaction of two 1/2 electron spins: 

"^KdaS-YrSar (17) 

occupying the overlapping orbitals (here 7 refers to that of the three d-orbitals 
which overlaps with the particular a-orbital) with the effective exchange con- 
stant given by: 

Kda = tl{/^E-Xa + ^E-Xd) > 0, (18) 

where 

/^Ed^a = ad-aa + Uaa " 2Udd + l^/AJdd = Aa + Uaa - 2Udd + 15/4Jdd > 0, 

AEa^d = aa-ad + SUdd + 50/4 J^^ = -Aa + Wdd + 50/4 J^^ > 0, 
Aa = ad - a„ > 0; AEa^d > AEd^a- 
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Aa - gad > 0, , 




= 4- 


Ad - gad > 0, 


Id 


- 7^° - 


Cd>0;l'd^a^d-{nd-l)Ud4 




- I'd- 


Udd 


la 


- <- 


■ C„ > 0. 



On the other hand the charge transfer energies can be expressed through the 
spectral ionization potentials of the respective ions, their electron affinities 
and the energy shifts Cd and of these quantities induced by the Coulomb 
field of the surrounding crystals already mentioned: 



(19) 



(here we omitted the intraatomic exchange parameters Jdd known to be by 
orders of magnitude smaller than other quantities relevant here and included 
the electron- hole interaction energies gad previously absorbed in a's). 

The interaction eq. (17) must be repeated for each interacting pair of 
electronic spins (pair of orbitals, coupled by the electron hopping operator). 
For the V ion in the crystal the terms appear for the ^- and ry-states on each 
atom. 

For the pair of metal ions interacting through the [TCNE]2^ unit the one 
electron hopping is effectively possible not only between the states in the t2g- 
manifolds of the ions involved, but also between other remaining (i-orbitals. 
By this the states admixed to the spin degenerate bare gound state of the 
pair of ions may have either lower or higher spin which leads both to the 
antiferromagnetic contribution of the form: 

^ (20) 
and of the ferromagnetic contribution of the form: 

2 

I Jdd (21) 



Kd 



\AEd^d, 

where the hopping parameters tdd and t'dd have the same order of magnitude. 
The terms of the antiferromagnetic sign appear for each pair of the d-orbitals 
coupled by the hopping, whereas the terms of the ferromegnetic sign appear 
all singly occupied orbitals in both shells provided an empty or a doubly 
occupied in the given rf-shell couples with a singly occupied orbital in other 
d-shell through the hopping. 
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B Phenomenological Hamiltonian for effec- 
tive spins and its relation to the spin Hamil- 
tonian 



The phenomenological Hamiltonian written in terms of the effective spins eq. 
(7) to be used for modeling the entire crystal has the form of eq. (8) 

In order to obtain it we follow the general recipes given in Ref. [16] for 
obtaining the spin-wave spectrum and write the general Heisenberg equations 
of motion for the spin-raising operators 5'^^, 'S'^'^. They read: 



^M^dr — J± {^dr^dr±h ^dT^dr±h) 



+ ■^11 ^dr {^ar + ^ar+a+c + ^ar+a + ^ar+c) 
~ ^dr ('S'ar + '^'ar+a+c + ^ar+a + '^'or+c) 
^Ir {^dr + ^dr-a-c + ^dr-a + '^'dr-c) 



z 

dr— a— c 



+ ^dr-a 



+ ^dr-c 



(22) 



In the spin-wave approximation the operators and are replaced by their 
average values in the magnetic (ordered) phase: 



2 ' ^ar ~^ (Si 



^dr ~^ {^d 

SO that the equations of motion get the form: 



1 

2 



''"'dt'-^ar 



-J± (fS", 



drib '-''-'dr 



+ J 



_i_ c+ -1-9+ 

ar+a+c ' ar+a ' ar+c 



3 /" C+ -1 C+ 
~\ [pdr + ^dr-a-c + ^dr-a + 



dr— c 



+ 2^+ 
- 65+ 



Going to the Fourier transforms of the raising operators we get: 



(23) 



(24) 



''"'dt'-^dk 
''"'dt'-'ak 



3 J_L (cos A';^ 
-J, 



2 '-'dk 



1) ^dk 

25+ 



^k5jk + 25t 



1 -|- exp(i/Ca) + cxp(i/cc) + exp(iA;a + ike 



(25) 



which must be complemented by analogous system of equations for the spin 
lowering operators 5^, 5^. 



24 



On the other hand the Heisenberg equations of motion for the operators 
S^j., 5'^, S^j., and S^^. as derived from eq. (6) (after the Fourier transformation 
is performed) form a system of four equations of motion for the Fourier 
components of the spin-| raising operators for each wave vector k. It can 
be rewritten with use of k-dcpcndcnt 4x4 matrices acting on the vectors 
with the components Sq^, 5'^, S^^^: 



where 



M(k)5+ 



(26) 



M(k) = A + T(k) + L(k) 



(27) 



and the matrix 



/ 4 



dd 



-2 

-2 



V 



-2 -2 ^^ 

4-2 

-2 4 

/ 



describes the spin fluctuations in the d-shells, the matrix 



(28) 



r(k) = K, 



dd 
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/ 1 
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) 



(29) 



describes the spin wave propagation in the 6-direction (transversal to the 
V-TCNE) planes and the matrix 



Lfk) 



K. 



da 






V 



with 



<5ck 





1 





1 

e' 





1 



+ e- 





2 Qck 

-2 / 



(30) 



(31) 
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desribes the spin- wave propagation in the ac-plane i.e. in the individual 
V-TCNE layer. Going to the linear combitations of the spin fluctuation 
operators in the d-shell: S^^ + + >S'^; —5'^ + S^^k, —Sq^ + 5*^ introduces 
the fluctuation of the effective spin of the (i-shell (the first combination) 
and incidentaly diagonalizes the sum of the first two matrix terms A + T(k) 
yielding the eigenvalues: 0; 3Km (cos k;, — 1) ; 6Jdd ~ ^K^d > of which the 
zero one refers to precession of the | spin in the acceptor orbital, the next one 
refers to precession of the total spin of the d-shell, and the doubly degenerate 
highest eigenvalue corresponds to excitations changing the total spin of the 
d-shell. The interaction matrix L(k) in this basis acquires the form: 



K. 



da 



I -12 

3Qok 
3Qck 
V 3 (g„k + <9ck) 



Q:k 

4 

-2 
2 



2Q, 



ak 



Qak 
-2 

4 
2 



2Qck QaV. Qck \ 



2 
2 

4 



/ 



(32) 



The states corresponding to the excitations of the d-shell (the above matrix 
is written so that they are the second and the third ones) are of much higher 
energy than the states corresponding to precession of the effective spins in two 
types of sites of the model. For that reason the former can be excluded. To do 
so we treat L(k.) as a perturbation with the smallness parameter and in 
the zero order we obtain an operator acting in the two-dimensional subspace 
spanned by the vector with the components Sq^ + -|- — S'^; 5'^'^: 



-^Kda 

h {Qak + Qck) 



K, 



da 



2I 



da 



-3Kdd{l- coskb)) 



(33) 



Comparing eq. (33) with eq. (25) and noticing that Qk = Qlk + Qck ^c arrive 



Si 



S: 



to the conclusion that they coincide after setting 
and Kdd = —J±', Kda = 3J|| which establishes the required re^ 
the spin wave treatments of the Hamiltonians eq. (6) and eq. 
order correction to eq. (33) has the form: 



_ 1 _ 1 I C2 I 

- 2 - 3 I'^dl 

at ion between 
(8). The first 



^da 



9J, 



dd 



I (2 — COS ka — COS kc + 2 cos kc cos ka) 

— J^k 




(34) 



but it is not used hereinafter. 
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C Mean Field treatment of the phenomeno- 
logical Hamiltonian. 

In order to obtain a structure dependent mean field estimate for the Neel 
temperature of the V-TCNE material we use the method described in Ref. 
[21]. Accordingly the Neel temperature for a system comprising two types 
of magnetic centers with effective spins Sd and Sa so that each center of 
the with spin Sd has Zdd neighbours of the same type, Zda neighbours with 
the spin Sa etc. with the interactions between the nearest neighbors of the 
specific type given by Jdd, Jda ^tc. satisfies the equation: 

JddZddSd {Sd + 1) — 20^^ JdaZdaSd {Sd +1) _ q 

JdaZadSa {Sa + 1) Jaa^aaSa {Sa + 1) — 

where we used the fact that J ad = Jda- For the structure represented in 
Fig. 1 and modelled by the Hamiltonian eq. (8) we set Sa = \- Then 
Zad — Zda — 4; Zdd — 2; Zaa — (acceptors do not have acceptor neighbors 
in this model). In the notation of eq. (8) Jdd — J±', Jda — J\\, so that we 
obtain: 

^MF ^ 2Jx ^ ^ ^■Iisj{Sd + lf + plSd{Sd + l) 

which in the limit of strong anisotropy yields the following mean field esti- 
mate: 

^MF ^ 2J^^^ + ;l) + J|| ^p,{Sd+l) (35) 

which fiows to the the limiting expression used in Ref. [8] for the single Fe- 
TCNE layer with the variance of the factor of two which results from the 
different definition of the Hamiltonian in Refs. [8,21]. 

D Spin-wave model of magnetic properties of 
V(TCNE)2 

Now we address the spin- wave theory of the M(TCNE)2 ferrimagnets de- 
scribed by the phenomenological Hamiltonian eq. (8) in order to derive its 
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magnetic properties. The temperature dependence of magnetization and the 
Neel temperature (that at each the magnetization of each sublattice vanishes) 
are controlled by the spectrum of the lowest energy excitations: spin waves 
(magnons). Calculation of these properties is customary performed with use 
of the Holstein-Primakoff Ref. [17] representation of the magnons. This latter 
had been many times derived for the ferrimagnets with relatively simple crys- 
tal lattice Refs. [18,19]. These derivations are based on the Green's function 
techniques. The present case differs from those described there by a com- 
bination of the alternation of the interaction sign (Jy and J± have opposite 
signs) and a relatively complex form of the structure factors which makes the 
Green's function too cumbersome and also prevents from using directly the 
the general formulae derived previously. We perform the required derivation 
for the present structure for certainty using the equation of motion method. 
It evolves as follows. Each term in the phenomenological Hamiltonian eq. 
(8) has the form: 



J12S1S2 — J] 



12 



(36) 



If the pair of the above spins couples ferromagnetically ( J12 < 0) the Holstein- 
Primakoff operators hi{h'l) annihilating (crearing) an elementary excitation 
(magnon) at the z-th site are introduced by the relations: 



Si = \j2Si-hibi 

S- = htpS,-n, (37) 
SI = Si-hi;hi^hjhi 

The operators hiihf) so defined obey the boson commutation relations: 



— Sii' 



(38) 



After being inserted in the above Hamiltonian term the square roots are 
expanded and the terms non higher of the second order in the boson operators 
bi are kept (linear spin-wave approximation) so that one gets: 



J12S1S2 = J12S1S2 - J12 {Sih2 + hiS2) + Ji2\J S1S2 (b^b2 + 6^61 



(39) 
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If the pair of effective spins couples antiferromagnetically (J12 > 0) the 
Holstein-Primakoff bosons are obtained by a more complex trick which, first, 
requires a transformation of the spins residing on the m-th site in the unit 
cell by the rotation matrices 00^' 



Si 



SiUJi] Si — uj^ Si 



(40) 



In the assumption of the common quantization axis for all spins in the crys- 
tallographic unit cell one can select the rotations as follows: 



UJi 












n 









1 









-1 













[0 








(41) 



in the basis of the coordinate spin components x,y,z or 





(1 









(0 


1 









1 


:) 


; (^2 = 


1 


















lo 








(42) 



in the basis of the tensor spin components +, — , z. Then the interaction term 
rewrites: 



J12S1S2 — —J12S1S2 + -J12 (^SiS2 + "S*! <S'2 ) 



(43) 



The boson operators for the transformed spins are then defined by the same 
relations as the nontransformed ones which yields the following interaction 
term: 



J12S1S2 = -J12S1S2 + J12 {Sifi2 + niS2) + Ji2\J S1S2 (btb2 + 6162) 



(44) 



where the anomalous products bfb'^ and 6162 appear. 

Described procedures apply to each pair of the interacting spins (Jy > 
0; Jj_ < 0) in the Hamiltonian eq. (8) yielding the following: 



HSW = Yl {-4:J\\SaSd + J±Sl + 4J|| {havSd + Safldr) - '2J±S, 



d^dr 



-\-J\\\j SaSd [-Dp {Ar + Ar+a + ^r+c + ^r+a+c) 



(45) 



A 



r+a 



A 



r+c 



A 



r+a+c 



+ J^Sd {D+Dr+^ + Z}++b^r 
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where D^{Dr) and A^{Ar) are the magnon creation (annihilation) operators 
at the respective sites in the crystal and the 'numbers of magnons' operators 
necessary to calculate the magnetization at each site are: 

har^ AtAr;har = D+Dr (46) 

Introducing the Fourier transforms of the magnon creation operators by the 
relations: 

Dt = ^Eexp(-ikr)D+; 



= -^Eexp(-ikr)A+ (47) 



and by the and the hermitean conjugate ones for the annihilation operators 
we obtain the spin-wave Hamiltonian in the form: 



+ J|l ^SaSd [fik^k^-k + ^^DiA\\ (48) 
+ 2 J^Sd {cos h-l)D+Dy,) 

where the unnecessary constant is omitted and the structure factors Qk co- 
incide with those defined by eq. (25). 

The above Hamiltonian produces the Heisenberg equations of motion 

for the annihilation operators and which are coupled with the analogous ones 
for the creation operators due to the presence of the anomalous terms: 



iMk - 4J\\SdA^ + J\\yfs^d^_kD: 



ihD+ = -4:JiiSaD+- Jii^SaSdniA_y, + 2 J^Sd {cosh- 1)D+ (49) 
ihi)^ = 4J\\SaD^ + J\\^/s^d^kA\-2 J^Sd {cosh -I) 



One can easily see that in the above system the first one is coupled only to 
the second whereas the third one is coupled only with the fourth. Thus eq. 
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(49) reduces to a pair of 2x2 matrix eigenvalue/eigenvector problems to be 
solved for the stationary magnons: 



4: Sd - g k ^/SaSd^l-k \ [ U]^^ 



(50) 



45'a — 2Sdak — £k V SaSdi^k ) [ ^k \ _ (-> /r-i \ 

^ -v^^k -^Sd-eJ[-yJ~ ^ ^ 



where 



Ok = (1 -cos/cb) (52) 

The matrix eigenvalue problems eqs. (50), (51) each yield two solutions for 
of which one is negative in either case and must be rejected (see Ref. [16]). 
The nonegative solutions of both problems represent the spectrum of magnon 
excitations having the form: 

4 = r k±{2{Sd-Sa) + Sdak) (53) 

Tk = \/{2{Sd + Sa) - Sda^f - SdSa \n^f 

where the " — " sign corresponds to the gapless band with the magnon anni- 
hilation operators G^, whereas the "+" sign in the expression for the energy 
corresponds to the excitations with a gap and the magnon annihilation op- 
erators Fk given by the relations: 

i^k = UkAy,-VkD\ (54) 
Gk = XkDi,-y^A\ (55) 

The site populations by the magnons are given by 

M = ^ E (^k^k) ; (udr) = ^ E (d^d^) (56) 

k k 



where 



2 







' + l?/kl 






+ l^kl 
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D+D^) = \vkr{FkF+) + \uknG+Gk) (57) 



and the (...) averaging is performed over the equihbium state of a ferrimagnet, 
which can be rewritten in comphance with the commutation rules for the 
boson operators: 



(^Gk)) 



-Dk-Dk) = l^kl^ (1 - {fiFu)) + \uk\^ (nck) 



The average population of a magnon state at a temperature T 
presses through the energy of the corresponding magnon 



^k 



J||£ 



(58) 
9/kB ex- 

(59) 



following the boson statistics: 



{flGk) 



1 -1 



exp 



e 



; {fiFk) 



exp 



k4 
9 



1 -1 



(60) 



In the low-temperature limit the contribution of the gap magnons is exponen- 
tially small so that {fipjs) can be set equal to zero. The expansion amplitudes 
UkjVkjXk, and yk are themselves immaterial and the only required densities 

I 2 I 2 

are |xk| and |yk| given by: 



(61) 



1 |2 ^ 
pkl 1 






bkl^ J 


Sd + Sa 





(here we neglected the k-dependence of the magnon amplitudes). 
The summation prescribed by eq. (56) replaces according to: 



-V — 



(62) 



BZ 



with the first Brillouin zone being a cube with the side of 27r. 

Evaluating the integrals is based on the advantage of the long wave ap- 
proximation for the energy. Two versions of the latter are relevant: in the 
first low-temperature regime <^ \ J±\ < J\\ the gaplcss branch of the magnon 
energy spectrum has the form: 



C2 
'^d 



Sd — Sa " Sd — Sa 



SdSg 

Sd — S, 



■Jiik 



(63) 
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which by the standard moves (see e.g. Ref [22]) brings us to the estimates 




(here C is the Riemann (^-function) . This is in some variance with Rcf. [18] 
since in the present case the structure factors for the magnetic interactions 
corresponding to the reahstic model of the material under study have been 
used rather the simple cubic ones and the definition of the exchange param- 
eters in the Hamiltonian differs by a factor of two (each pair of interacting 
spins counts once here so that the coefficient of the interaction parameter in 
the diagonal matrix element of the equation of motion equals to the num- 
ber of neighbours of the correcponding type). Taking into account quadratic 

5 

corrections to the densities eq. (61) yields the higher order corrections oc 9^. 

3 

Neglecting these corrections one gets the Bloch T2 law eq. 9 for the sponta- 
neous magnetization with the critical (Neel) temperature (one at which the 

spontaneous magnetization disappears and which in the present model coin- 
cides with the point where the sublattice magnetizations disappear either) 
to be found from the equiation: 

1 {S,i-Sa)-^ /3\ 



C(^)— ^ = l• 



Expanding the magnon energy upto the second order in all components 
of the wave vector is possible only if the temperature is the smallest enery 
scale. In the case of very high anisotropy the second low temperature regime 
is possible when |J_l| -C ^ -C Jy. In this case the gapless branch of the 
magnon energy spectrum has the form: 

= (/^J||^^-2-^^Jx(l-coskb)) (65) 



1 2 1 2 

and the integration of I Xk I or j^kj divided by 



first performs 

(following Ref. [20]) over the planar projection kp of the wave vector k, which 



exp ( ^ 1 - 1 
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reduces to 



ab 



/ 



47r2y [exp(pA;2)-a 



(66) 



where 



exp 



Sd — Sa 9 



[1 — coskb) 



SdSa J\\ , 

I P — -] = 

Sd — Sa Sd + 



After substituting z — the integration yields: 



log(l - a) 



(68) 



Provided a is an exponential function with a small negative argument, log(l- 
a) replaces by the logarithm of the absolute value of the argument: 



log(l — a) log 



1^1 - cos kb) , 



^Sd — Sa 9 

so that the expression to be intrgrated over kb becomes 



1 Sd — Sn 



StT^ SdSa J\\Sd + Sa iSa 



S,i 



log 



|Jx 

Sd — Sa 



(1 — COS kb) 



which immediately yields 

1 Sd — Sa 



1 



Sd 



A'K SdSa JwSd + SaiSa 

The sublattice magnetizations are: 

I Sd — Sa 1 



log 



q2 



Sd — Sa 9 



Sd 

Sn 



1 



47r Sd + J| 



^ 1 (Sd — Sa 

log 



C2 



Jl 



(69) 



(70) 



(71) 



so that the Neel temperature is obtained by the condition of the evanescence 
of the sublattice magnetizations which in the present approximation coincide 
with that of the total spontateous magnetization. 
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Table 1: Saturation magnetization Mg at 2 K; ordering temperature and 

effective exchange energy between the metal residing effective spins 

in the assumption of the mean field connection between Tc and for 

M(TCNE)2 after [6]. 



Metal 


Saturation 


magnetization Mg 


Tc in K 


Jeff in K 




emu Oe mol-1 


Spins per M atom 






ya 


10 300 


2 


~ 400 


53 




6 100 


1 


^ 400 




Mn 


19 000 


4 


107 


6.1 


Fe 


16 900 


3 


121 


10 


Co 


8 000 


~ 1.5 


44 


5.9 


Ni 


15 800 


3 


44 


11 



« The HM form. 
The LM form. 
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Figure 1: Structure of Fe(TCNE)2 as coming from the synchrotron radiation 
study [8]. Each unit cell contains four formula units. The solvent molecules 
are omitted for clarity. 
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Figure 2: Hypothetical structures of V(TCNE)2 following [9] (left) and following [10] (right). 
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Figure 3: Schematic representation of the intermediate structures between 
the quadrupled "principal" (left upper corner) and the experimental (right 
lower corner) ones. 
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Figure 4: V-Densities of states in two spin channels for the hmiting and an intermediate structures 
between the quadrupled "principal" (left) and the experimental (right) ones. Observe different scale of 
the abscissa for different pictures. The magnetization in units of number of unpaired electrons per unit 
cell is indicated on top of each picture. 
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Figure 5: N-Densities of states in two spin channels for the hmiting and an intermediate structures 
between the quadrupled "principal" (left) and the experimental (right) ones. Red/green is the DoS of 
the N atoms in the V-TCNE layers, blue/yellow is that of the "to be dangling" ones; olive green/dark 
red is that of those involved in the V-V bonding through the [TCNEjg" units. Observe different scale of 
the abscissa for different pictures. The magnetization in units of number of unpaired electrons per unit 
cell is indicated on top of each picture. 
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Figure 6: C-Densities of states in two spin channels for the hmiting and an intermediate structures 
between the quadrupled "principal" (left) and the experimental (right) ones. Red/green is the DoS of 
the C atoms in the V-TCNE layers, blue/yellow is that of those forming the C-C bonds in the [TCNE]!^ 
units. Observe different scale of the abscissa for different pictures. The magnetization in units of number 
of unpaired electrons per unit cell is indicated on top of each picture. 
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Figure 7: The projection of the DoS in two spin channels to the [TCNE]2~ 
units, in the "experimental" structure [8]. Red/green is the DoS of the N 
atoms bonding the V-TCNE layers; blue/yellow is that of the "dangling" 
ones; olive green/dark red is that of the C atoms. 
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